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Abstract. The magnetic insulator Gd2Sn207 is one of many geometrically frustrated 
magnetic materials known to exhibit a nonzero muon spin polarization relaxation 
rate, A(T), down to the lowest temperature (T) studied. Such behaviour is typically 
interpreted as a significant level of persisting spin dynamics (PSD) of the host material. 
In the case of Gd2Sn 2 07, such PSD comes as a surprise since magnetic specific heat 
measurements suggest conventional gapped magnons, which would naively lead to 
an exponentially vanishing A(T) as T — > 0. In contrast to most materials that 
display PSD, the ordered phase of Gd2Sn207 is well characterized and both the 
nature and the magnitude of the interactions have been inferred from the magnetic 
structure and the temperature dependence of the magnetic specific heat. Based on this 
understanding, the temperature dependence of the muon spin polarization relaxation 
through the scattering of spin waves (magnons) is calculated. The result explicitly 
shows that, despite the unusual extensive number of weakly dispersive (gapped) 
excitations characterizing Gd2Sn 2 0j, a remnant of the zero modes of the parent 
frustrated pyrochlore Heisenberg antiferromagnet, the temperature dependence of the 
calculated A(T) differs dramatically from the experimental one. Indeed, the calculation 
conforms to the naive expectation of an exponential collapse of A(T) at temperatures 
below ~ 0.7 K. This result, for the first time, illustrates crisply and quantitatively the 
paradox that presents itself with the pervasive occurrence of PSD in highly frustrated 
magnetic systems as evinced by muon spin relaxation measurements. 

1. Introduction 

Muon spin rotation (/xSR) [U [2] is a local probe that is sensitive to the magnetic 
environment at the microscopic level in condensed matter. Briefly, fiSR experiments 
involve injecting positive muons polarized along some direction into a condensed matter 
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sample. Each muon comes to rest within the sample and its spin precesses in the local 
magnetic environment. It may also be subject to processes that cause it to depolarize. 
The muon is an unstable particle and decays on a time scale of the order of //s. Since 
it is mediated by weak interactions, the decay cross section in the muon rest frame 
is anisotropic in space with the anisotropy axis depending on the muon polarization 
direction at the time of the decay. One therefore measures the asymmetry of the decay 
with respect to the initial polarization axis to obtain information about the effect of the 
material on the magnetic moment of the muon. 

One of the many uses of /zSR is in probing static and dynamic processes in magnetic 
materials [H [2] . Consider, for example, a conventional insulating magnetic material that 
undergoes a transition to long range order at some temperature T c . For temperature 
T > T c and in zero applied magnetic field, the material is a paramagnet with no static 
moments. Injected muons therefore experience only fluctuating fields which bring about 
an exponential relaxation of the asymmetry. In the critical region near T c , the muon 
relaxation rate increases, diverging at T c as the spins exhibit critical slowing down. 
Below T c , muons precess in the local magnetic field. If the magnetic environment is 
homogeneous, the asymmetry exhibits oscillations at the Larmor frequency of the muon. 
However, samples are often inhomogeneous and there is also generally a distribution of 
muon stopping sites, so the oscillations decay or are altogether undetected because the 
asymmetry probes the precessional dynamics averaged over the whole distribution. Such 
relaxation due to static fields, leaves a residual asymmetry which can decay via genuine 
dynamical processes in the sample. In conventional magnets, with long range order or 
with frozen moments, as in spin glasses, this dynamical relaxation rate A should tend 
to zero as the temperature vanishes. 

In contrast, /xSR studies have uncovered, over the last two decades, many materials 
in which there is anomalous relaxation rate at low temperatures that is temperature 
independent. Among them are quasi two dimensional [31 EJ El El E] and 3D magnets 
[9]- [25] . Persistent spin dynamics (PSD) has been detected both in long range ordered 
magnetic systems P [101 HH H21 [131 HH H3 US HI] , in candidate spin liquids [El HHJ ED] 
and systems with freezing transitions with no apparent long range order [201 EH 1221 [23] . 
Also, there is persistent dynamics in the low temperature spin ice state of the Ising 
magnet Dy 2 Ti 2 0y although the temperature dependence has not yet been explored 
[21 [25]. 

To date, there is no quantitative microscopic theory that accounts for PSD in 
any material in which the phenomenon is seen. The problem is especially paradoxical 
in magnets with long range order because spin dynamics should be frozen out as the 
temperature is lowered well below the phase transition temperature. In this article, we 
highlight the problem by focussing on a geometrically frustrated magnetic material 
for which the magnetic structure has been characterized, the principal interactions 
identified and for which the low energy excitations can be computed reliably. The 
material in question is the rare earth pyrochlore Gd2Sn 2 0y [T5], [161 ESI E3 EEl [291 130] . 
another material that exhibits PSD [151 EE 121]. In the next section, we introduce this 
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material in some more detail with an outline of its phenomenology and a description 
of the available microscopic theory for Gd2Sn 2 0y. Section [3] is a discussion of the field 
distribution in Gd2Sn2C>7 and an estimate of the muon position and distribution of 
fields that the muons experience. In Section IH we take the spectrum of excitations 
and compute from it the relaxation rate one would expect below T c and compare with 
the experimental rate, finding significant discrepancy. As far as we are aware, this is 
the first detailed microscopic calculation of the muon spin depolarization rate, A, in a 
material exhibiting PSD. Finally, we present in Section |5] a comparison of our results 
with the experimental data and close with a general discussion of PSD in magnetic rare 
earth pyrochlore oxides. 

2. Gd 2 Sn 2 7 

Gd2Sn207 is a magnetic insulator with space group Fd3m in which the Gd 3+ ions 
occupy the sites of a pyrochlore lattice. The magnetic levels have S = 7/2 and L = 
so the crystal electric field has an effect only through the mixing of states outside these 
L = levels. The material has a negative Curie-Weiss temperature #cw = —8.6 K 
and undergoes a first order transition at about 1 K to long range order [30J with a 
Palmer- Chalker magnetic structure [31] and an ordered moment of about 7/xb [30J . 
This magnetic structure strongly suggests that the dominant magnetic interactions 
are antiferromagnetic exchange with additional (perturbative) magnetostatic dipolar 
interactions. These observations as well as, for example, specific heat [26] and electron 
spin resonance (ESR) measurements [28] have allowed for a determination of the 
exchange coupling as well as a more refined picture of the microscopies to include a weak 
XY-like single ion anisotropy. The leading contribution to the single ion Hamiltonian is 

if aniso = | A| ^ ' ^af (1) 

where \x labels the face-centred cubic space lattice sites and a runs over the tetrahedral 
primitive basis of the pyrochlore lattice. The local sublattice z a axis is oriented in the 
local [111] direction on sublattice a. The anisotropy coupling A has been estimated to 
be about 0.14 K [32]. Further contributions to the crystal field of Gd2Sn207 are smaller 
and have been neglected [32]. The magnetic interaction Hamiltonian is 

Hint = Jcx'y^ J Jj • Jj +^ >r n n y^[J» " Jj ~ 3(Jj • Rjj)(Jj ■ Rjj)]|Rjj| 3 (2) 

where is the position vector of ion % with respect to ion j. The estimate of J ex = 0.27 
K which comes directly from 6qw [26] provides an adequate fit to the ESR data in 
Ref. [28] which measures the uniform contribution to the spin wave spectrum. The 
dipolar coupling V = (/i /47r)(g/z B ) 2 /r^ n = 0.049 K where r nn = \/2a/4 is the distance 
between nearest neighbor Gd 3+ ions and a is the unit cell edge length of 10.46A [26J and 
g — 2 is the Gd 3+ Lande factor. One can also incorporate exchange couplings beyond 
nearest neighbour. An estimate for these couplings exists [33] but we neglect them here 
because they do not significantly affect our results. 
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The first experiment to explore the electronic spin dynamics in Gc^S^C^ was a 
Mossbauer spectroscopy experiment which concluded that the Gd 3+ moments continue 
to fluctuate at least down to 27 mK. The two /zSR experiments of Bonville et al. [15] 
and Chapuis et al. [16J report asymmetry data below T c with two clear oscillation signals 
for time t < 0.3/^s corresponding to a pair of local fields of 207(1) mT and 442(1) mT. 
The initial measured asymmetry below T c in Ref. [16] is reported to be only about one 
quarter of that above T c which is perhaps due to muons at other stopping sites and 
experiencing larger internal fields giving precession signals which are outside the (time) 
resolution of the spectrometer. The relaxation rate A of the nonoscillating part of the 
asymmetry shows, in both experiments, as the temperature is lowered, a clear and sharp 
peak at T c followed by a plateau at about 0.6/is -1 extending down to the lowest observed 
temperature of 0.02 K. 

3. Field Distribution in Gd 2 Sn 2 7 

The internal fields determined from the //SR spectrum in two experiments [T5J [16] 
are very similar suggesting that the muon location is rather well-defined in Gd2Sn 2 07. 
Both experiments find a satisfactory fit to the asymmetry oscillations using two distinct 
fields. This is perhaps unexpected since, even allowing for a fixed muon location, there 
are three distinct Palmer- Chalker magnetic domains (up to time reversal). In this 
section, we explore the magnetic fields in the vicinity of the oxygen ions, where the 
positive muon is expected to hydrogen-bond, assuming that (i) the only field that the 
muons experience comes from the static Gd 3+ magnetic moments and that (ii) the muon 
does not produce a local lattice distortion. Taking a cubic cell of volume (2a) 3 with 
magnetic moments in the Palmer- Chalker state [31], we compute the magnetic field 
within spherical shells of inner radius I A and outer radius 1.5 A about each oxygen ion. 
We find locations where the three possible symmetry-related Palmer-Chalker magnetic 
domains give rise to two different fields within experimental precision which have the 
same magnitudes observed experimentally. Specifically, in units of a, one of these 
locations is (0.4402, 0.5005, 0.5627) where the three Palmer-Chalker domains produce 
fields of 206 mT, 207 mT and 441 mT. Since two of the fields are identical to within 
experimental precision, one would expect that the weight that field contributes to the 
oscillatory part of the asymmetry would be twice that of the other. However, the analysis 
of the asymmetry [16J finds equal weights for the two oscillatory components. 

4. Calculation of A 

In this section, we compute the relaxation rate of the muon due to a two-magnon Raman 
process. This is identical to calculations of spin lattice relaxation in NMR [31]. The 
muon down spin (|J,)) to up spin (|f)) transition rate is given by the Fermi Golden rule 

W = -fr E<K f ;t |^ t -Gd|U)| 2 )th^f -w!-7m#) ( 3 ) 

a,fe,i,f 
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where a sum is taken over all initial (i) and final (f) magnon states and ^^H is the muon 
Zeeman splitting in the local field H. (. . .} t h denotes a thermal average. 

The reference (unperturbed) Hamiltonian of Gd2Sn2C>7 is the sum of Eqns. (CQ) and 
02]) and the interaction between the magnetically ordered moments and the muon spin 
I is an interaction between dipole moments given by 

iVod = (£) 9^Bl,hJ2 jjjL {(I ■ *S a (R")) - 3 (i ■ R£) (5S a (Rn -k)}. 

fi,a 

We have included only the deviation of the moments from their ordered (Palmer- 
Chalker) directions 5S because only the fluctuation part can lead to dynamical relaxation 
of an individual muon. The static part of the interaction Hamiltonian, which is not 
shown, causes the muon to precess in the local field and a static distribution of such 
internal fields produces the nondynamical relaxation described in the Introduction. To 
proceed, we have to consider the possible processes that might give rise to relaxation. We 
work in the linear spin wave approximation using Holstein-Primakoff bosons [281 1331 135] . 
The simplest process is one where a single magnon is created or annihilated, accompanied 
by a muon spin flip. Such a process appears to lowest order from the transverse spin 
components, for these are linear in the magnon creation and annihilation operators but 
it is forbidden by energy conservation because the spin wave gap is much larger than the 
muon Zeeman splitting per unit field. To illustrate this point, we compute the Zeeman 
splitting of a muon to be hj^ = 0.0065 KT _1 . Since the local field is of the order of 10 2 
niT [151 HE], the Zeeman energy is about 6.5 x 10~ 4 K, while the spin wave gap is about 
0.5 K in Gd 2 Sn 2 7 [281 EH |35], so a single magnon typically has an energy that exceeds 
the muon spin flip energy by three orders of magnitude. For this reason, we neglect 
the transverse Gd 3+ spin fluctuations and drop the Zeeman energy from the energy 
conserving delta function in Eqn. ([3]). With only longitudinal fluctuations contributing 
to the relaxation, let 

r(Rn = (£) sj^fcl { (I ■ Zfl ) - 3 (i • R£) (z fl • } (4) 

where unit vectors z a are along the local quantization axis of the Gd 3+ spins and depend 
only on a because the ordering wavevector is q = in the Palmer- Chalker state [31 J . 
The interaction can now be written as 

iVGd=-x>»( RM )n R/1 ) ( 5 ) 

l-i.a 

with h a = cj l (R At )c a (R /1 ) for = S — 5S* = S — c\c a in terms of Holstein-Primakoff 
boson operators. The lowest order energy conserving contribution to the relaxation rate 
is a two magnon Raman process whereby the muon spin is flipped via the scattering of 
a magnon. The matrix element for this process is 

M aaW (k,k') = (rv(k) + l,n 6 /(k') - 1;| \H^ Gd \n a (k),n b (k');i) (6) 

where n a is the boson occupation number for magnon branch a, with a running over the 
four magnetic sublattices, a = 1, ... ,4, of the Palmer- Chalker state [31 J . 
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Figure 1. Panel 1(a) shows the spin waves frequencies, cj a (k) for GcbS^Cv along 
paths in the Brillouin zone of the BCC reciprocal lattice of the FCC space lattice shown 



in Panel 1(b) The magnon gap never reaches its minimum of about 0.7 K along this 
path. 



The calculation of the spin wave spectrum for this material has been discussed 
elsewhere [26] 12H1 E3j [35] and we refer to the literature for further details. We replace 
the spins in the Hamiltonian H = if an i SO + H int with Holstein-Primakoff bosons and 
expand in powers of 1/S up to and including powers of S° which gives a harmonic 
Hamiltonian. One finds that the terms linear in the boson creation and annihilation 
operators, c* and c, vanish in the Palmer- Chalker spin configuration, as required of a 
semi-classical ground state of the Hamiltonian H. The c\c operators are transformed 
via a canonical transformation Q ab (k), an 8 x 8 matrix, into spin wave operators and 
a. The transformation, 

4 

c a (k) = ^gfa 6 (k) + Qf +4 aS(-k) 

6=1 

4(-k) = J2Ql +4b a b (k) + Q^al(-k), 

6=1 

where index a runs over the sublattices, is such that the spin wave operators also satisfy 
Bose commutation relations and transform the spin Hamiltonian into the form 

H = H (o) + 5 £ E Wa (k) + E ^( k ) K(kK(k)) , 

k a k a 

where the spin wave frequencies are u a (k). The spectrum along line segments between 
points of high symmetry in the Brillouin zone is shown in Fig. [1] where the canonical 
transformation has been computed numerically [33], [35] . 
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We can now compute the matrix element .M ao /&6/(k, k') of Eqn. (jBJ) using 
a+(k)|n(k)> = ^(k) + l|n(k) + 1) and a(k)|n(k)) = y^kj\n(k) - 1). We find 

.M aa w(k,k') = 5 aa ,5 w [(Qi) ca Q$ + (Q\>) cb+4 Q-T} 

fl,C 

x v/n a (k) + fv^(k0exp(2(k - k') • R£)<t |/ C (R")| I) 

where Q* ab is the complex conjugate of Q ab . Computing the squared modulus of this 
matrix element, we make the approximation that the phase factor exp(z(k — k') • (R£ — 
R2)) averages to zero unless R^ — R^ = so that 

l^l 2 ~ $> b (k')K(k) + i)F c (q*»q£ + Q*t+ 4 Q-+ 4 ) {Ql a Qi> b + Q-W 4 ) , 

c 

where F c = \{f |/ C (R M )| l)\ 2 . The transition rate W in Eqn. ([3]) then becomes 

W =lTc^[ dk I rfk'Vn 6 (kO(n a (k) + l)% a (k)-o; b (kO)F4(g-gr)(Q c k 'Qk' 6 ) 
ft "bz ^bz ./bz a 6 c 

+(Q-W4 a+4 )(Qlk' 4 e 6 ^ 4 ) + (Qk ca Qlk +4 ) (Qk'Q-k' 4 ) + (Q-^Q?)(Q-V~ 4 Qv)} , 

where the thermal average at inverse temperature /3 of the boson occupation numbers 
is the distribution n a (k) = (exp(/3a; a (k)) + l) -1 . The integrals can be regrouped into a 
single energy integral by introducing densities of states, 

<(u;) = 1 jf dkQ k 6 QJf iV 2 c6 (u,) = 1 jf rfkQi b + 4 Q^ +4 

JV?( W ) = l jf dkQM 4 iVf ( W ) = 1 J dkQ^Q? = JV*», (7) 

where the u indicates that the integrals are evaluated for fixed energy and the indices b 
and c run over sublattices 1 to 4. The normalization factor M is chosen to ensure that 
(1/jV) J dk = 1. Finally, the transition rate W is given by 



^ ' m=l a,6,c 



(8) 



It remains for us to evaluate F c from Eqn. (TjJ. Because the Gd2Sn 2 07 sample 
in Refs. [151 HH] is a powder, we must average over the crystal orientations with 
respect to the initial muon polarization axis. For this purpose, it is convenient to 
define three separate coordinate systems: (i) the muon polarization system fit,, (ii) 
the spin polarization system on each sublattice S a and (iii) a system about the axis 
R^. We introduce rotation matrices and V^f such that (S a ) a = fiip and 
(S a ) a = V°P(Tk%)p. These allow to write / C (R^) as 

/C(R " } = I (£) 9 ^ B ^ h w? { {u? " lUr) " 3V; " ^ v % ( u < x " ^ 

Performing the angular average, we obtain 

FC = 5 (£) 2 £ |R^ + 3 ^ ( (F - )2 + W + 
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Figure 2. Panel (a) is the two magnon Raman muon A(T) for GcbS^Or and Panel 
(b) is a plot of the densities of states (DOS) for each of the four branches of the spin 
wave spectrum. The lowest branch in the density of states has a low energy tail that 
extends down to the spectral gap of about 0.7 K. 



We evaluate Eq. (JSJ) averaged over the three possible Palmer-Chalker magnetic domains 
and over all crystal orientations. The densities of states N^(u) in Eqn. (j7|) were 
computed by Monte Carlo integration by introducing energy bins over the range of 
magnon energies. For points k distributed uniformly in the first Brillouin zone, 0J a (k) 
are determined whereupon the integrand in the density of states is added into the 
corresponding energy bins. The muon location was taken to be the position given in 
Section [3j Finally, the muon polarization relaxation rate A = 2W, where the factor 
2 arises from the fact that both spin-up and spin-down populations contribute from 
the relaxation rate [M], is shown in Fig. [2] plotted as a function of temperature up 
to the transition temperature. Compared to the experimental plateau at 0.6/is^ 1 , 
the two magnon Raman relaxation rate is very small - about one-hundredth of the 
experimental relaxation rate in the vicinity of T c . In addition, the relaxation rate is 
strongly temperature dependent. A calculation of the relaxation rate for a gapped 
cubic lattice antiferromagnet [31] gives a trend A oc AoTexp(— A/T) which is illustrated 
for our problem in the inset to Panel (a) of Fig. [2j The T — » limit of A(T) in the inset 
is consistent with the spectral gap A ~ 0.7 K. 

5. Discussion and Conclusions 

We have computed the relaxation rate of the muon spin polarization due to the scattering 
of a pair of magnons in the ordered phase of Gd 2 Sn 2 7 . As far as we are aware, this 
is the first such detailed microscopic calculation of /iSR A(T) for a frustrated magnet 
exhibiting PSD. We confirm explicitly that the rate is strongly temperature dependent, 
in line with the observation that the magnon density should fall off exponentially at 
temperatures less than the spectral gap A w 0.7 K. The two-magnon relaxation rate 
is a strongly temperature dependent background to the experimentally measured A. 
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However, in spite of the strong temperature dependence, we find that A just below T c 
is smaller than the experimental rate by more than two orders of magnitude. Since 
the variation of the computed rate with temperature is smaller than the experimental 
resolution, the two-magnon (Raman) process in Gd 2 Sn 2 7 is completely hidden by 
whatever mechanism is responsible for the A ~ 0.6 /is -1 relaxation rate plateau. 

The relaxation of the muon polarization of the order of 0.6 /zs -1 is therefore entirely 
inexplicable over the whole temperature range below T c « 1 K within the usual context 
of spin lattice relaxation. This result puts aside the possibility that the cause of this 
phenomenon comes from unusual features of the spin wave spectra of geometrically 
frustrated magnets [TU [36]. Higher order magnon processes are also suppressed at 
temperatures below the spin wave gap (37] so neither the plateau nor the A ~ 0.6 /zs -1 
rate can be explained by considering the scattering of conventional (gapped) magnons. 
It is interesting, in light of this result, to recall details of some suggestive experimental 
findings. 

Firstly, PSD is typically regarded as a distinctive feature of magnetic materials 
on frustrated lattices as opposed to conventional low temperature magnets. However, 
beyond this feature, there are not many similarities between the materials exhibiting 
anomalous low temperature relaxation. Notably, one finds PSD in materials without 
long range order (e.g. Tb 2 Ti 2 07 [HI [19] and Yb 2 Ti 2 7 [20]). Indeed, PSD uncovered by 
/iSR is often taken to be a good evidence from a local probe of spin liquid behaviour. It is 
not clear that the mechanism underlying the spin dynamics is the same across the large 
spectrum of materials exhibiting PSD. However, if the mechanism is common to all these 
materials, then its manifestation in collective paramagnetic or spin liquid phases may 
constitute further evidence that magnons are irrelevant to low temperature A plateaux. 
An interesting finding is that the A plateau is little affected by doping Tb 2 Ti 2 7 with 
nonmagnetic impurities [19], suggesting that the plateau, at least in this compound, is 
not due to conventional collective excitations. To emphasize, the observation of PSD in 
(gapless) spin liquid systems is not our immediate concern. Rather, it is the occurence 
of PSD in pyrochlores such as Tb 2 Sn 2 7 [121 H31 HH HS1 EH1 EHl SQl SI], Er 2 Ti 2 7 
[9], Gd 2 Ti 2 7 pH QI] Gd 2 Sn 2 7 [151 EE] as well as spinels [UJ such as FeAl 2 4 and 
CoAl 2 04 that display long range order as exposed by neutron scattering experiments 
that we believe is particularly paradoxical. 

One might wonder whether the muon itself might be contributing to this 
phenomenon. In this context it is worth commenting on Tb 2 Sn 2 7 . In this material, 
spin dynamics below T c has been observed in /^SR [121 UHl EES] as well as with a neutron 
probe [HHl [39l [40], [41] , albeit on much shorter timescales from those in /zSR experiments 
and apparently with a significant temperature dependence [10]. In conclusion, the 
microscopic origin of the ubiquitous phenomenon of PSD in highly frustrated magnetic 
materials remains an important and unsolved problem. 
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